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Q Abstract 

Though extensively studied, hardness, defined as the resistance of a material to deformation, still remains a challenging issue for a 
formal theoretical description due to its inherent mechanical complexity. The widely applied Teter's empirical correlation between 
1^ hardness and shear modulus has been considered to be not always valid for a large variety of materials. Here, inspired by the 
classical work on Pugh's modulus ratio, we develop a theoretical model which establishes a robust correlation between hardness 
and elasticity for a wide class of materials, including bulk metallic glasses, with results in very good agreement with experiment. 
The simplified form of our model also provides an unambiguous theoretical evidence for Teter's empirical correlation. 
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1. Introduction 

, Despite the great efforts, to understand the theory of hard- 
ness and to design new ultrahard materials is still very chal- 
lenging for materials scientistslfll 0, Q |4|]. During the past few 
years, several semi-empirical theoretical models llH lil I?! S 0] 
have been developed to estimate hardness of materials based 
on: (0 the bond length, charge density, and ionicity |5], (//) the 
strength of the chemical bonds |6], (///) the thermodynamical 
concept of energy density per chemical bonding [7], and (iv) 
the connection between the bond electron-holding energy and 
hardness through electronegativity [8], and (v) the temperature- 
dependent constraint theory for hardness of multicomponent 
bulk metallic glasses (BMGs) |9]. Experimentally, hardness is 
a highly complex property since the applied stress may be de- 
pendent on the crystallographic orientations, the loading forces 
and the size of the indenters. In addition, hardness is also char- 
acterized by the ability to resist to both elastic and irreversible 
plastic deformations and can be affected significantly by defects 
(i.e., dislocations) and grain sizes [10]. Therefore, hardness is 
jiot a quantity that can be easily determined in a well-defined 
absolute scale It has been often argued [13] that hardness 
measurements unavoidably suffer of an error of about 10%. All 
these aspects add huge complexity to a formal theoretical defi- 
- nition of hardness lUQ, H H . @1 ■ 

Within this context, to find a simple way to estimate hard- 
ness of real materials is highly desirable. Unlike hardness, the 
elastic properties of materials can be measured and calculated 
in a highly accurate manner Therefore, it has been histori- 
cally natural to seek a correlation between hardness and elas- 
ticity. The early linear correlation between the hardness and 
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bulk modulus (B) for several covalent crystals (diamond. Si, Ge, 
GaSb, InSb) was successfully established by Gilman and Cohen 
since 1950s [10, 11]. Nevertheless, successive studies demon- 
strated that an uniformed linear correlation between hardness 
and bulk modulus does not really hold for a wide variety of 
materials ifTii [H [Tsll . as illustrated in Fig. [TJa)- Subsequently, 
TeterjT^j established a better linear correlation between hard- 
ness and shear modulus (G), as illustrated in Fig. [Ub). This 
correlation suggests that the shear modulus, the resistance to 
reversible deformation under shear strain, can correctly pro- 
vide an assessment of hardness for some materials. However, 
this correlation is not always successful, as discussed in Refs. 
ifsi flsifTill . For instance, tungsten carbide (WC) has a very large 
bulk modulus (439 GPa) and shear modulus (282 GPa) but its 
hardness is only 30 GPa ifisll . clearly violating the Teter's lin- 
ear correlation [see Fig. [lib)] [Q. Although the link between 
hardness and elastic shear modulus can be arguable, it is certain 
to say that the Teter's correlation grasped the key. 

In this manuscript, following the spirit of Teter's empirical 
correlation, we successfully established a theoretical model on 
the hardness of materials through the introduction of the classic 
Pugh modulus ratio of G/B proposed in 1954 [T^. We found 
that the intrinsic correlation between hardness and elasticity of 
materials correctly predicts Vicker's hardness for a wide variety 
of crystalline materials as well as BMGs. Our results suggest 
that, if a material is intrinsically brittle (such as BMGs that fail 
in the elastic regime), its Vicker's hardness linearly correlates 
with the shear modulus (Hy - 0.151 G). This correlation also 
provides a robust theoretical evidence for the famous empirical 
correlation observed by Teter in 1998. On the other hand, our 
results demonstrate that the hardness of crystalline materials 
can be correlated with the product of the squared Pugh's mod- 
ulus ratio and the shear modulus (//,, = 2(A;^G)"^^^ - 3 where 
k is Pugh's modulus ratio). This formula provides the firm evi- 
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Figure 1: (Color online) Correlation of experimental Vickers hardness (//,,) 
with (a) bulk modulus (B) and with (b) shear modulus (G) for 39 compounds 
(Tables |2] and [3). Inset of panel (b): Hy vs. G for 37 BMGs (see Table [T). 
The solid line denotes empirical Teter's fitting values, whereas dashed lines 
correspond to the value derived from Eq. (|6). The black and hollow squares 
denote data taken from Refs. ( |12t il||). 



dence that the hardness not only correlates with shear modulus 
as observed by Teter, but also with bulk modulus as observed 
by Oilman et al. Our work combines those aspects that were 
previously argued strongly, and, most importantly, is capable 
to correctly predict the hardness of all compounds included in 
Teter's ||l2t]. Oilman's flUSl, Gao's Jj and Simunek's 16] sets. 
Also, our model clearly demonstrates that the hardness of bulk 
metallic glasses is intrinsically based on the same fundamental 
theory as the crystalline materials. We believe that our relation 
represents a step forward for the understanding and predictabil- 
ity of hardness. 

2. Model and Results 

According to Vicker ifioll . the hardness of //,, is the ratio be- 
tween the load force appUed to the indenter, F, and the indenta- 
tion surface area: 



2Fsin(0/2) 
^2 ' 



(1) 



where d and 6 are the mean indentation diagonal and angle be- 
tween opposite faces of the diamond squared pyramid indenter, 
respectively (Fig. |2|i. In order to derive our model, we first 
assume that (;) the diamond squared pyramid indenter can be 
divided into four triangular based pyramid indenters and that 




Figure 2: (Color online) Illustration of indentation in terms of the squared di- 
amond pyramid indenter. The red framework highlights one of four triangular 
based pyramid indenters. 



(;/) the Vicker's hardness is measured within the elastic scale. 
Then, for each triangular based pyramid, one can define the 
shear modulus G as. 



4Atan(a) 



(2) 



which specifies the ratio between shear stress and the shear 
stain. In terms of our model the exact shear area A on which 
the shear force (F) acts is unknown. But, the deformation area 
A* [A*-^cfitan(a)] delimited by the klO' triangle is well de- 
fined by the indentation geometry. Therefore, we can express 
the exact shear area (A) as: 



A = cA* = —d tan(a) 



(3) 



where c is the proportional coefficient. It is clear that under 
elastic shear deformation the deformation area (A*) will be ex- 
tremely small. However, upon real hardness measurements the 
deformation area (A*) should be large enough so that the coeffi- 
cient c can be safely neglected and A ^ A*. Under this assump- 
tion, the equation (|2| can be revised as following. 



2F 



(4) 



d^tan^{a) 

Combining equations ^ and (|4|, the Vicker's hardness reads 

H, = Gtan^{a)sin(ei2) = 0.92Gtar?(a) (5) 

where the term sin(6'/2) is intrinsically determined by the in- 
denter itself, which can be considered as a constant (originated 
from the Vicker's hardness, see equation ([!}). For the diamond 
squared pyramid indenter with - 136° then sin(0/2) is equal 
to 0.92 for Vicker's hardness measurement. In an ideal form of 
indentation, tan(a') - 0.404 because of a = (c.f.. Fig. |2]l. 
Therefore, the equation (|5]l can be simplified as, 

//, = 0.151G. (6) 



2 



Equation (|6]l represents a robust theoretical evidence of 
the linear correlation behavior observed by Teter [|l2i . as re- 
flected by the data shown in Fig. [Ub)- Residual discrep- 
ancies should be mainly attributed to the neglection of plas- 
tic deformation effects. Remarkably, we found that Eq. (|6]l 
is also valid for BMGs. Using the experimental shear modu- 
lus G = 38.6 (36.6) GPa |17] for Pd4oNi4oP20 BMG, the es- 
timated Vicker's hardness is 5.83 (5.53) GPa, in consistency 
with the experimental value of 5.38 GPa |17]. Similarly, for 
Fe4iCo7Cri5Moi4Ci5B6Y2 BMG by using the experimental G 
(84.3 GPa lUi) we obtained = 12.7 GPa, in nice accord 
with the measured Vicker's hardness of 12.57 + 0.22 GPa 1 18]. 
As illustrated in the inset of Fig. [Ub), the agreement with the 
experimental values is highly satisfactory for all 37 BMGs col- 
lected here (see Table [TJ- Considering that BMGs are brittle 
materials without plastic fractures the data in the inset of Fig. 
[lib) strongly convey that our proposed formula (Eq. |6) is in- 
trinsically connected to the shear modulus of materials if they 
fail in an elastic regime. 
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Figure 3: (Color online) Experimental Vicker's hardness as a function of the 
product k^G (k=GIE). All data ai'e collected from literature (see Tables |2]and 



It is highly difficult to realistically take plastic deformation 
into account in our modeling scheme. However, the indentation 
after a real hardness measurement shows the permanent plastic 
effect, which is, of course, reflected by the ratio [namely, 
equal to tan(Q')] (see Fig. |2]i. Note that the depth of the in- 
dentation, loo', is parallel to the direction of shear deformation. 
We reasonably assume that its size should be closely correlated 
with the shear modulus of G, whereas the expansion wideness 
of the indentation, Ipo", is perpendicular to the direction of the 
loading force, hence, with almost little connection to the shear 
deformation. Therefore, the expansion wideness seems to re- 
flect the ability to resist to compression effects, a property that 
should be related to bulk modulus, B. Accordingly, we pro- 
posed the following relation. 



Finally, combining the equations (|5]) and (|7]i, the hardness can 
be written as. 



oc G(G/By 



(8) 



It is interesting to note that in Eq. (jSJ the ratio of G/B is the fa- 
mous modulus ratio proposed by Pugh in as early as 1954 lll6ll . 
In his pioneer work, Pugh derived that the strain at fracture can 
be measured as e oc (B/G)^. Indeed, hardness can be defined as 
the resistance to the applied stress at the critical strain of e that 
the system can sustain before yielding it to fracture. This clearly 
provides fundamental support for our model (Eq. |8). Impor- 
tantly, Pugh also highlighted a relation between the elastic and 
plastic properties of pure polycrystalline metals and stated that 
G/B is closely correlated to the brittle and ductility of materials 
UM- the higher the value of G/B is, the more brittle the mate- 
rials would be fl^. Otherwise, the materials are expected to 
deform in a ductile manner with a low G/B value. This relation 
has been extensively accepted and a ppli ed not only to metal but 
also to high-strength materials ^25, 26, 27|. In principles, the 
covalent materials (such as diamond and c-BN) have the highest 
hardness but they are obviously brittle with a larger Pugh mod- 
ulus ratio. The strong covalent bonds indeed create a significant 
resistance to initialize the plastic flow to pin the dislocation, re- 
sulting in a quite high hardness. Conversely, ductile materials 
with a low Pugh's modulus ratio are characterized by metallic 
bonding and low hardness. It is thus highly reasonable to estab- 
lish a correlation between hardness and the modulus ratio G/B 
in Eq. (O. Thus, we revise further the Eq. (|8) in the following 
form. 



H, = cr'G"; {k = G/B) 



(9) 



where //,,, G and B are the hardness (GPa), shear modulus 
(GPa) and bulk modulus (GPa), respectively. The parameter k is 
the Pugh's modulus ratio, namely, k = G/B. C is a proportional 
coefficient. In order to derive the parameters C, m and n, we 
first selected ten materials with diamond-like (diamond, c-BN, 
j6-SiC, Si and Ge), zinc -blende (ZrC and AIN) and rock-salt 
structures (GaP, InSb and AlSb) because their hardness, bulk 
and shear moduli are well-known (see Tables |2]and|3]l. By an- 
alyzing these data we found that C=1.887, m=1.171 and n - 
0.591. Hence, Eq. (|9]) becomes. 



H, 



1.887fc''^iG°-5^' 



lM7(k^G)°-^^\ 



(10) 



In order to assess the validity of this relation, we plot in Fig. 
[3]i/,, against k^G for a series of compounds. These data show a 
clear and systematic trend with k^G and firmly establish a direct 
relation between hardness and k^G. By fitting the data of Fig. [3] 
and revising further Eq. (fTOt we arrive to the final formula: 



2(k^G)'''' 



(11) 



tan(a) oc G/B 



(7) 



from which we see that the hardness correlates not only with 
the shear modulus but also with the bulk modulus. Physically, 
the bulk modulus only measures isotropic resistance to vol- 
ume change under hydrostatic strain, whereas shear modulus 
responses to resistance to anisotropic shear strain. Although the 
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Table 1: The comparison between the predicted Vicker's hardness by Eq. (6) [Hy = 0. 151G] and the experimental data for 37 bulk metallic glasses are shown, 
together with experimental Young modulus (E) and Shear modulus (G). The [XX];[YY] in the first column denotes the reference numbers: XX is the reference for 
elastic constants and YY is the reference for the Vicker's hardness. 
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bulk modulus was thought to be less directly connected with 
hardness ||5|, the Pugh's modulus ratio k clearly contributes 
to the Vicker's hardness. Equation (flU demonstrates that, if 
bulk modulus increases, hardness would decrease as long as the 
shear modulus remains unchanged, and vice versa. This behav- 
ior can be understood by the fact that if the Pugh's modulus ra- 
tio G/B gets smaller with increasing bulk modulus, the material 
would become more ductile. Its hardness can be thus expected 
to have a lower value. Taking the example, the compounds 
TiN and jS-SiC have almost the same experimental shear moduli 
lHIH of 187.2 GPa and 191.4 GPa, respectively. However, 
the experimental bulk modulus of TiN (318.3 GPa)||2i ^ is 
larger by about 42% than that of y6-SiC (224.7 GPa). In terms 
of our formula, y6-SiC is found to be harder than TiN, in agree- 
ment with the experimental observationsjsl |3^ [Expt (Calc in 
this work): TiN with H, = 23 (20) GPa and /3-SiC with H, = 34 
(33) GPa]. 




10 20 30 40 50 60 70 80 90 100 
Experimental H (GPa) 



Figure 4: (Color online) Congelation between experimental and theoretical 
Vickers's hardness (//,,) for 39 compounds, as compared with the estimated 
data from the models llHl (see Tables|2]and|3). 

To further assess the performance of our model (Eq. [TT] we 
show in Fig. |4] a comparison between the estimated and ex- 
perimental values for a series of compounds (see Tables |2] and 
O, confirming a good agreement. Also WC, which is wrongly 
found to be a superhard (49 GPa) material within Teter's lin- 
ear correlation, is now predicted to have a Vicker's hardness 
of 29.3 GPa in very good accordance with experimental value 
(30 GPa fl3]). In paiticulai-. Figs. |3]and|4] convey that our 
proposed formula reproduced very well the Vicker's hardness 



bU6LBC2N 



0331 IMI). 



for all well-known superhard materials (Daimond 
iBBIllll, c-BN i,!!, C-BC5 f32^, and y-Bzs 

The interesting case is the compound of ReB2, which was 
thought to be superhard Il82ll. Although its Vicker's hardness 
was debated extensively 1821 l83l 18411 . there is now a wide- 
accepted consensus that its Vickers' hardness of 30. 1 GPa at 
the large loading force of 4.9 N 11451 14611 . Using the experi- 
mentally measured bulk and shear moduli (43\ (B - 273 GPa 
and G = 382 GPa) and in terms of our Eq. (flU . the Vicker's 
hardness is derived to be 32.9 GPa, in nice agreement with the 
experimental data l45ll . 



Another attention has to be paid to the case of BeO. Using the 
experimental bulk and shear moduli {B = 230 GPa and G = 206 
GPa) [85], its Vicker's hardness is calculated to be 36.7 GPa 
within our current model. This value is well within the scale of 
the experimentally measured results from 32 to 38 GPa ll42l 85 j 



for polycrystalline boron suboxide sintered samples, although a 
Vicker's hardness of 45 GPa was reported for the single crystals 
under a loading force of 0.98 N ll42ll . Indeed, the light loading 
force of 0.98 N is not large enough to obtain a real hardness. 
It is thus expected to have a smaller hardness if a loading force 
larger than 0.98 N is applied. Our estimated value for BgO 
is also in good agreement with the derived value of 37.3 GPa 
through a very recent thermodynamic model of hardness [86.1 . 

We further estimated two more phases of carbon (C^ and M- 
carbon), which were suggested to be superhard 
91]. Utilizing elastic shear and bulk moduli obtained in Ref. 
1 89], the Vicker's hardness of C4 is calculated to be 69.0 GPa 
(c.f.. Fig. O that is comparable to the superhard c-BN. More- 
over, through using the calculated bulk and shear moduli (B = 
415 GPa and G = 468 GPa IH) for the M-carbon phase, we 
obtained its Vicker's hardness of 8 1 .0 GPa (c.f.. Fig. |3]l, placing 
M-carbon in between BC2N and diamond, agreein g w ell with 
the value (83.1 GPa) obtained by Simunek's model 08811 . 

In addition, from Fig. |4]all estimated data are in good agree- 
ment with those obtained from pervious models [5, 6]. Never- 
theless, we would like to emphasize that, although our proposed 
model can reproduce well the results obtained by Cao's f3\ 
and Simunek's [6] models, the underlying mechanism is sub- 
stantially different. Gao's and Simunek's models are based on 
bond properties such as bond-length, charge density, ionicity 
and their strengths and coordinations in crystalline materials. 
Differently, our model depends totally on the so-called poly- 
crystalline moduli (bulk and shear modulus as well as Pugh's 
modulus ratio), which indeed response directly to the abilities 
of resistance under loading forces for polycrystalline materials. 
As demonstrated above, for polycrytalline materials the intro- 
duced Pugh's modulus ratio in our model plays a crucial role in 
elucidating plastic deformation, which is intrinsically different 
from all known semi-empirical hardness models dH S HI . 

3. Discussion and Remarks 

The hardness of a material is the intrinsic resistance to defor- 
mation when a force is applied [1]. Currently, a formal theoreti- 
cal definition of hardness is still a challenge for materials scien- 
tists. The need for alternative superhard and ultrahard materials 
for modern technology has brought a surge of interest on mod- 
eling and predicting the hardness of real materials. In particular, 
in recent years several different semi-empirical models for hard- 
ness of polycrystalline covalent and ionic materials have been 
proposed. Gao's model is mainly based on bond length, charge 
density, and ionicity [5]. Simunek's model employs the strength 
of the chemical bonds and its framework in crystalline materials 
[6]. Mukhanov's model utilizes the thermodynamical concept 
of energy density per chemical bond [7]. Li's model is mainly 
based on the bond electron-holding energy hardness through 
electronegativity |@]. Smedskjaer's model correctly predicts 
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Table 2: Comparison between theoretical values within our current model and experimental values as compared with available theoretical findings obtained through 
the models of Gao 1 5] and of Simunek 1 6]. Furthermore, the bulk and shear moduH are compiled in this table. In the last column, "e" and "c" denote elastic data (G 
and B) from direct experimental and theoretical calculations, respectively. The Pugh's modulus ratio k is compiled in this table. For details, see text. The [XX];[YY] 
in the first column denotes the reference numbers: XX is the reference for elastic constants and YY is the reference for the Vicker's hardness. 
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iinwiek 




Uiamond |z8J;|5J 


CO c c 

535.5 


442.3 


1 oil 

1.21 1 


ACT 

95.7 


96 


AO C 

93.6 


AC A 

95.4 


e 


Uiamond |z8J;|5J 


C /I 

548.3 


A C 

465.5 


1 1 TO 

1.178 


AO A 

93.9 


96 


AO C 

93.6 


AC A 

95.4 


c 


Diamond |3jJ;|jJ 


COA 


/I 1 A 

431.V 


1 AC 

1.205 


AO C 

93.5 


C\£. 

96 


AO C 

93.6 


AC A 

95.4 


c 


Diamond 1 12] 


CO C A 

535.0 


443.0 


1 AO 

1.208 


AC A 

95.4 


A/C 

96 


AO C 

93.6 


AC A 

95.4 


e 


BC2N [36];[37J 


446. U 


A AT A 

403.0 


1 1 AT 

1.107 


T/I A 

76.9 


T/C 

76 


T 

78 


T 1 A 

71.9 


c 


BC2N [12];[37J 


/I /I C A 

445.0 


A AO A 

408.0 


1 AA 1 

1.091 


T C A 

75.4 


T/C 

76 


T 

78 


T 1 A 

71.9 


e 


BC5 [38];[32] 


'~iC\ A A 

394.0 


376.0 


1 A /I 

1.048 


£.£ T 

00.1 


T 1 

71 






c 


T> XT r 1 . rci 

c-BN [lo;!^ 


/I AC A 

405.4 


/I AA A 

400.0 


1 A 1 /I 

1.014 


c 

65.2 


66 


64.5 


ci 

63.2 


e 


c-BN [28J;|5J 


A AT /I 

403.4 


A AT T 

403.7 


A AAA 

0.999 


/CO 

63.8 


66 


C A C 

64.5 


CI 

63.2 


c 


„ "D XT noi.rci 

c-bN |z6];l^ 


3oz.z 


3 /j. / 


1 Ai n 
1.01 / 


/;o 1 
63.1 


CC 

66 


CA C 

64.5 


CI 

63.2 


c 


„ 13 XT r A AT . rci 
C-BN |4U];|j] 


A(\A A 

404.4 


/I A 

384.0 


1 ACO 

1.053 


£.0 

68.2 


CC 

66 


64.5 


CI 

63.2 


c 


T) XT r 1 . rci 

c-BN fl2];[5J 


A AA A 

409.0 


A AA A 

400.0 


1 AOO 

1.023 


66.2 


CC 

66 


C A C 

64.5 


/CO 

63.2 


e 


y-B28 f41];[33J 


236.0 


224.0 


1 AC /I 

1.054 


Ac\ r\ 

49.0 


C A 

50 






c 


BeO [12];[42J 


A/1 A 

204.0 


228.0 


A AC 

0.895 


A 

36.4 




38 






e 


yS-SiC f28];[5] 


191.4 


224.7 


0.852 


32.8 


34 


30.3 


31.1 


e 


y8-SiC [28J;q] 


196. 


/I A 

224.9 


A OT A 

0.874 


1 A C 

34.5 


/I 

34 


A 

30.3 


01 1 

31.1 


c 


y8-SiC [43J;|5J 


1 AA 

190.2 


AA 

209.2 


A AAA 

0.909 


C C 

35.5 


/I 

34 


A 

30.3 


01 1 

31.1 


c 


y8-SiC |44J;[5J 


1 O/C c 

186.5 


220.3 


0.846 


A 

32.0 


/I 

34 


A 

30.3 


01 1 

31.1 


e 


y8-SiC [12];[5J 


1 A/C A 

196.0 


226.0 


A CI 

0.867 


1 A 1 

34.1 


/I 

34 


A 

30.3 


01 1 

31.1 


e 


S1O2 [12];|12J 


'^OA A 

220. 


T AC A 

305.0 


A TO 1 

0.721 


OA A 

29.0 




33 


OA A 

30.4 




e 


ReB2 [45J;|46] 


273.0 


TOO A 

382.0 


A T 1 C 

0.715 


A 

32.9 


OA 1 

30.1 






e 


WC [47J;|5J 


OA 1 

301.8 


A10 C\ 

438.9 


f\ £.00 

0.688 


00 A 

33.4 


OA 

30 


26.4 


01 C 

21.5 


e 


WC [9;[5J 


252.0 


/I OA A 

439.0 


A /I 

0.642 


OA 

29.3 


OA 

30 


26.4 


01 C 

21.5 


e 


B4C [48];[12] 


1 AO A 

192.0 


226.0 


A OCA 

0.850 


000 

32.8 


A(3 

30_ 






e 


VC [Inis workj;!^ 


AA 1 

209.1 


AC C 

305.5 


A £.0 C 

0.685 


O/C 

26.2 


OA 

29 




OT 

27.2 


c 


ZrC |49|;|50| 


169.7 


223.1 


A n £ 1 

0.761 


O/I 

26.3 


c 

25.8 






e 


ZrC |51|;|50J 


1 O'^ C 

182.5 


228.3 


A TAA 

0.799 


OA A 

29 A 


c 

25.8 






c 


ZrC [51J;[50J 


1 C A 

185.9 


A 

228.0 


A 1 C 

0.815 


A C 

30.5 


c 

25.8 






c 


ZrC [52j;|5a] 


1 £.C\ C 

169.6 


223.3 


A TCA 

0.759 


O/I 

26.2 


c 

25.8 






e 


ZrC 1 12|;|50] 


1 /I/^ A 


223.0 


r\ n A A 

0.744 


c 

25.2 


c 

25.8 






e 


liC |4yj;|oJ 


lo2.2 


/I A 


A TCI 

0. /53 


OT 1 

2 /.I 


/I T 

24. / 




100 
18.8 


e 


liC |j3J;|oJ 


1 T/:: A 

1 /o.y 


OCA 

250.3 


A nf\n 
0. /O/ 


/I c 

24.5 


/I T 

24. / 




100 
18.8 


c 


liC p4J;|oJ 


1 AO 


00/; A 

286.0 


A £C\'~1 

0.693 


oc 
25.8 


/I T 

24. / 




100 
18.8 


c 


nr-\ rcci.r/ci 

TiC [55J;|^ 


1 OT 

187.8 


0-11 T 

241.7 


A TTT 

0.171 


000 
28.8 


/I T 

24.7 




100 
18.8 


e 


iiC U2j;|oj 


1 00 A 

loo.O 


/I 1 A 

241.0 


A TOA 

0. /80 


OA A 

29.0 


/I T 

24. / 




100 
18.8 


e 


liN |5oJ;|50J 


183.2 


000 A 

282.0 


0.650 


00 C 

22.5 


00 
23 




1 T 

18.7 


c 


TiN [57J;[50J 


187.2 


T 1 T 

318.3 


A COO 

0.588 


OA A 

20.0 


00 

23 




1 T 

18.7 


e 


TiN f58J;[50J 


AC 

205.8 


OA /I 

294.6 


A £.C\C\ 

0.699 


O/C T 

26.7 


00 

23 




1 T 

18.7 


c 


TiM r^Qi-r^ni 
111% |jyj,|jL^ 




JZ.O.J 


U.OJ) / 


Zj.O 


Aj 




lo. / 


c 


RUO2 [60];[5J 


142.2 


251.3 


0.566 


15.7 


20 


20.6 




c 


RUO2 [61J;i:5J 


173.0 


248.0 


0.698 


23.7 


20 


20.6 




c 


AI2O3 


161.0 


240.0 


0.671 


21.5 


20 


20.6 




c 


AI2O3 


160.0 


259.0 


0.618 


19.2 


20 


20.6 




c 


AI2O3 [62J;|;5J 


164.0 


254.0 


0.646 


20.7 


20 


20.6 




e 


AI2O3 U2];[a 


162.0 


246.0 


0.659 


21.1 


20 


20.6 




e 


NbC \^\\^ 


171.0 


333.0 


0.513 


15.5 


18 




18.3 


c 


NbC[52];[&] 


171.7 


340.0 


0.505 


15.2 


18 




18.3 


e 


AIN [40]; [5] 


134.7 


206.0 


0.654 


18.4 


18 


21.7 


17.6 


c 


AIN [6^;L5J 


130.2 


212.1 


0.614 


16.5 


18 


21.7 


17.6 


c 



''B4C was suggested to be very hard in Ref. 1 79]. Mukhanov et al recently predicted that the Vicker's hardness of B4C was 45.0 GPa 1 86] in agreement with 
the reported experimental data of 45 GPa in Ref. ]80] (see Table 1 in Ref. |86]). However, we also noted that Teter IT2I1 ever summarized the hardness of B4C 
with a value of 30±2 GPa. In addition, the experimental value of 32-35 GPa was recently summarized in the handbook 1*81]. Therefore, here we quoted the 
experimental Vicker's hardness of 30 GPa, as summarized by Teter in Ref. 1 121 . 
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Table 3: (Table[2]continued) 



V^VJiiiUVJ U-llLld 


rj 




1- 






ncao 


S imunek 






1 7^ ^ 

i Z, J . J 


707 5 


5Q4 


1 5 7 


1 8 
i o 


91 7 


17 6 
1 / .u 


c 




132.0 


211.1 


0.625 


17 


18 


21.7 


17.6 


g 


AIN n 71-1^1 


1 78 


70^ 


6^1 
W.O J i 


1 6 Q 


1 8 
1 


91 7 

z, i . / 


1 7 6 
1 / .o 


6 


iNDiN lO/ I.IOI 


1 SS Q 
i JJ.7 


7Q7 


S^4 
O. J Jt- 


1 ^ 4 


1 7 




IQ ^ 
1 7. J 


6 


NhN n ?i-rfii 


1 56 


315 


495 


1 3 9 

1. J .y 


17 




1 9 5 

i.y 




HfN r78i-rsni 




315 5 


591 

\j .^y 1 


70 


17 






f 


HfN r281-r501 


164.8 


278.7 


0.591 


18.4 


17 






Q 




1 OS 7 


1 7S q 
i / J.J/ 


SQ8 
W. J70 


1 7 
1 J. / 


1 S 1 

1 J. 1 


1 8 1 
10.1 


1 8 S 

1 O.J 


6 


riaN n 91-rsi 


1 70 


710 


S71 

W.J / 1 


14 1 


1 S 1 

1 J. 1 


1 8 1 
10.1 


1 8 S 
i o. J 


C 


7rn^ rzim-rsi 


oo.w 


1 87 


471 


8 4 


1 T. 


1 8 
IW.O 




c 


z^ji w 2 1 ^ ^ J 1 1 J 


93.0 


187.0 


0.497 


9.5 


13 


10.8 






i31 IH? I.I J 1 


66 6 
uu.u 


Q7 Q 


680 


1 1 8 
i 1 . o 


1 9 


1 ^ 6 
1 .u 


1 1 ^ 
1 1 . J 


c 


Si r7ni-rsi 
1 /ui.i J 1 


64 


Q7 Q 
y 1 .y 


654 


10 Q 


1 9 


1 ^ 6 


1 1 ^ 
1 1 . J 


c 


"?i r7m-rsi 


6^ 7 


QO 7 


6Q7 


1 1 8 


1 9 


1 ^ 6 
1 J.o 


I 1 ^ 

I I .J 


c 


Si r7n-rsi 


61 7 
O i . / 


Q6 ^ 
70. J 


640 


1 7 


1 9 


1 ^ 6 
1 J.o 


I 1 ^ 

I I .J 


c 


Si r7ii-rsi 


61 7 


89 


693 


1 1 5 

1 1 .J 


12 


1 3 6 

1 J .U 


1 1 ^ 
1 1 . J 




GaP rs2-r'si 


55.7 


88.2 


0.631 


9.3 


9.5 


8.9 


8.7 


g 


>jdr 1^;^ I.I J 1 


55 8 


88 8 


678 




Q 5 

y 


8 Q 

O . " 


8 7 

o . / 




yjdr 1^/ I.I Jl 


S6 1 

JO. i 


88 6 
00. o 


6^^ 

W.O J J 


q 4 


".J 


8 q 

o.y 


8 7 
0. / 




GaP r711-r'Sl 


61 9 


89 7 


690 


1 1 5 

1 1 ..J 


9 5 

y 


8 9 


8 7 

O . / 


^ 


AlP r7zll-lSl 
^ii 1 /^l.l Jl 




ou. w 


n S7n 

W. J / w 


7 1 
/ . 1 


Q A 


Q 6 

".U 


7 Q 


c 


AlP r7'^i-rsi 

/\lr 1 / Jl.l Jl 


SI 8 

J i .0 




n S7S 

W. J / J 


7 S 

/ .J 


Q A 


Q 6 

".U 


7 Q 


c 


AlP r7Si-rsi 

rtlr 1 / J I.I J 1 


4.8 8 


OU. V/ 


S67 

V. JU / 


7 
/ .w 


Q A 


Q 6 

".U 


7 Q 


c 


lIliN lOO I.IOI 




1 7^ Q 


AAA 


S 1 
J . 1 


Q 


10 A 


8 7 
o .z 


c 


TnN r741-rfil 

ILirH 1 /H-I.IUI 


77 




557 


Q 7 
y . 1 


Q 

y 


10 4 


8 7 
o .z. 


g 


Gp r70i-rfii 


53 1 


72.2 


736 


1 1 3 


8 8 


1 1.7 


9 7 

y . 1 


g 


Gp r7fllTfil 


43 8 


60 3 


726 


9 5 


8 8 


1 1.7 


9 7 


g 


rinAs r47i-rsi 

»J<lf\» IHI l.l J 1 


T-U — J 


7S 


670 
u.uzw 


7 8 
/ .o 


7 S 
/ .J 


8 


7 A 




ria A« r79i-rsi 




7S S 
/J.J 


610 

v/.U 1 " 


7 8 
/ .o 


7 S 
/ .J 


8 


7 A 




Ha As 14Q1-1S1 


Af^ 7 


7S A 


61 Q 


7 8 

/ .0 


7 S 
/ .J 


8 
o.w 


7 A 


c 


V^O, 1781-lSl 


77 S 

/Z.J 




A^7 


6 ^ 

U.J 


7 S 
/ .J 


7 7 




c 


V^O, 1781-lSl 


f^7 7 
uz. / 


IH-O. J 


A78 


J.J 


7 S 
/ .J 


7 7 




c 


Y^n, 17fil-lSl 


UU. J 


1 AQ ^ 


AAS 


6 
u.u 


7 S 
/ .J 


7 7 




6 


TnP r77i-rsi 




71 1 
/l.l 


A8^ 


^ 8 
J . o 


S A 


6 


S 1 

J . 1 


c 


TnP r47i-rsi 




77 S 

/Z.J 


A7S 

VJ.H- / J 


^ 6 

J.U 


S A 


6 
o.u 


S 1 
J. i 


6 


AiAs r77i-r'^i 


AA 8 


77 Q 

1 1 .y 


S7S 

U.J / J 


6 7 
u. / 


J 


8 S 
O.J 


6 8 
U.o 


6 


AiAs r79i-rsi 


AA 


78 ^ 

/ O. J 


S6Q 

U. JU" 


6 S 

U.J 


J 


8 S 
O.J 


6 8 
u.o 


6 


GaSh r771-rfil 






607 

v/.UU / 


S 8 


4 S 


6 

U.U 


S 6 


g 


riaSh 14Ql-lfil 
UdOD 1 ^? 1 . 1 O 1 


^A 1 


Sf^ A 


606 
w.uwu 


S 8 
J.o 


A S 


6 

U.U 


S 6 

J.U 


C 


riaSVi 1471-161 
UdOD 1^/ I.IOI 


^A ^ 


JU. J 


608 


S 8 
J.o 


A S 


6 
u.u 


S 6 

J.U 


C 


AlSh 178MS1 




S6 1 


S61 

v/.«JU 1 


4 6 


4 


4 Q 


4 Q 


g 


r\10L) 1 /ZJ,| JJ 




S8 9 
Jo.z 


SAQ 


A S 




A Q 


A Q 


c 




^7 S 

JZ. J 


J". J 


SA8 

U. Jt-O 


A 6 


A 


A Q 


A Q 




AlSb 1471-151 


31.9 


58.2 


0.549 


4.5 


4 


4.9 


4.9 


g 


TnAs [721-151 


29.5 


57.9 


0.509 


3.6 


3.8 


5.7 


4.5 


g 


InAs [49j;[5J 


29.5 


59.1 


0.499 


3.4 


3.8 


5.7 


4.5 


g 


InSb [47]; [5] 


23.0 


46.9 


0.490 


2.4 


2.2 


4.3 


3.6 


e 


InSb [72J;[5J 


22.9 


46.5 


0.492 


2.4 


2.2 


4.3 


3.6 


e 


InSb ]49};\5} 


22.9 


46.0 


0.498 


2.5 


2.2 


4.3 


3.6 


e 


ZnS [49]; [6] 


32.8 


78.4 


0.418 


2.5 


1.8 




2.7 


e 


ZnS |47J;[6J 


31.5 


77.1 


0.408 


2.3 


1.8 




2.7 


e 


ZnSe [47J;[6J 


28.8 


63.1 


0.456 


2.7 


1.4 




2.6 


e 


ZnTe [471; [6] 


23.4 


51.0 


0.459 


2.1 


1 




2.3 


e 


ZnTe \49}-\6} 


23.4 


51.0 


0.459 


2.1 


1 




2.3 


e 
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the hardness of multicomponent BMGs through temperature- 
dependent constraint theory_||9| . It has been noted that all the 
above models S 0, 0, H Ml have two major limitations: (i) 
Each one provides a satisfactory description only for a spe- 
cific type of materials: Smedskjaer's model treats BMGs with 
high degree of structural and topological disorder, whereas all 
other methods are essentially applicable to crystalline materials 
only, and (ii) they depends on different theoretical assumptions. 
Since these methods can only be used to predict the hardness 
of some specific materials a unified and general theory capable 
to account for the hardness of any material is still missing. If 
we look back the history of hardness, it can be easily found 
in literature that many scientists (Oilman, Cohen, Pugh and 
Teter) ESIIll Il6l II2II have tried to create a correlation be- 
tween the hardness and elasticity (a well defined quantity) since 
1950s. The most successful empirical correlation was proposed 
by Teter in 1998 11211 . who suggested that the hardness shows 
a quasi-linear correlation with the shear modulus. However, all 
these empirical correlations between hardness and bulk modu- 
lus (or shear modulus) turned out to be not fully successful. The 
main reason is that hardness indeed is a characteristic of a per- 
manent plastic deformation, whereas the elasticity corresponds 
to the reversible elastic deformation. Therefore, there seems to 
be a general consensus on the fact that the hardness in general 
does not depends neither on the bulk modulus nor on the shear 
modulus for polycrystal materials. Indeed, these correlations 
were heavily debated in recent years. 

In 1954 Pugh has proposed a relation between the elastic and 
plastic properties of pure polycrystalline metals and stated that 
the Pugh's modulus ratio (k - G/B) represents a good criterion 
to identify the brittleness and ductility of materials 1J_6J . It was 
found that material with a large k behaves in a more brittle man- 
ner, and that the higher the value of k is, the more brittle the 
materials are [,16il . On the other side, materials with a low k 
are expected to deform in a more ductile way. Basically, as ev- 
idenced from Tables|2]and|3]we found that the Pugh's modulus 
ratio, k, can be correlated with hardness. The hardest material, 
diamond, has the highest k value of about 1.2 and all widely 
accepted superhard materials have a highly large k value larger 
than 1.0. In addition, from Tables|2]and[3]one can see that with 
the progressive decrease of the Vicker's hardness the k values 
get progressively smaller Unlike the moduli of G and B, which 
only measures the elastic response, the Pugh's modulus ratio 
seems to correlate much more reliably with hardness because it 
responses to both elasticity and plasticity, which are the most 
intrinsic features of hardness. 

Therefore, through the introduction of the classic Pugh mod- 
ulus ratio proposed in 1954 fl^ and following the spirit of 
Teter's empirical correlation, we have constructed a theoreti- 
cal model of hardness. We proposed a new formula to calculate 
Vicker's hardness, - 2(fe^G)"^*^ - 3, for polycrytalline ma- 
terials. The most important aspect of our formula is that it cor- 
rectly predict the hardness of all compounds dataset considered 
in several recent models ||1 El H S E just under the 
condition of knowing the corresponding bulk and shear mod- 
uli. Furthermore, we proposed that, if a material is intrinsically 
brittle (such as BMOs), its Vicker's hardness linearly correlates 



with the shear modulus (//,, = 0.151 G). On the one hand, this 
simplified form provides a robust theoretical evidence why the 
famous empirical quasi-linear correlation observed by Teter in 
1998 are right for some materials. On the other hand, we found 
for the first time the Vicker's hardness can be linearly corre- 
lated with the shear modulus for BMGs (see Table [T]l. This is 
somehow unexpected because for BMGs there exists a univer- 
sal correlation between the Yough's modulus and the Vicker's 
hardness as documented in Refs. 

Finally, we still want to point out that our model (Eq. [TTJ 
may be not accurate to predict the hardness of pure metals (or 
metallic -bonding dominated materials with a highly low Pugh's 
modulus ratio). For instance, the hardness of fee Al is estimated 
to be 1.3 GPa, which is significantly larger than the experimen- 
tal Vicker's hardness of 0.167 GPa. This is mainly due to the 
fact that ductile metals can locally accumulate plastic defor- 
mation prior to fracture, which has not been considered in our 
model. 

In summary, via the aid of Pugh's modulus ratio, our work 
provides the firm evidence that the hardness not only correlates 
with shear modulus as observed by Teter (l^, but also with 
bulk modulus as observed by Oilman and Cohen 
By retaining the fundamental aspects of the previous proposed 
models, our model clearly demonstrates that the hardness of 
BMGs is intrinsically based on the same fundamental basis 
as the crystalline materials. Given the fact the elastic bulk 
and shear moduli can be accurately calculated by the state-of- 
the-art first-principles calculations, we believe that our finding 
is important for the community to design and develop ultra- 
hard/superhard materials and high-performance high-strength 
structural materials. 
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